The Hubble rate in averaged cosmology 
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The calculation of the averaged Hubble expansion rate in an averaged perturbed Friedmann- 
Lemaitre-Robertson- Walker cosmology leads to small corrections to the background value of the 
expansion rate, which could be important for measuring the Hubble constant from local observations. 
It also predicts an intrinsic variance associated with the finite scale of any measurement of Ho, the 
Hubble rate today. Both the mean Hubble rate and its variance depend on both the definition 
of the Hubble rate and the spatial surface on which the average is performed. We quantitatively 
study different definitions of the averaged Hubble rate encountered in the literature by consistently 
calculating the backreaction effect at second order in perturbation theory, and compare the results. 
We employ for the first time a recently developed gauge-invariant definition of an averaged scalar. 
We also discuss the variance of the Hubble rate for the different definitions. 



I. INTRODUCTION 

The late time Universe is not perfectly homogeneous 
and isotropic on all scales. The overdensities and voids 
that develop via gravitational collapse make it signif- 
icantly inhomogeneous. As a result, the notion of a 
maximally symmetric background geometry, which is the 
very basic foundation of the standard concordance model 
needs to be taken with extreme caution. Specifically, 
one would like to construct an averaged model which 
can suitably describe the Universe on sufficiently large 
scales, as a coarse-grained version of the actual distribu- 
tion of matter and energy in the Universe. In the last 
decade, this issue has attracted considerable attention in 
cosmology, in particular through the so-called averaging 
problem (see [1] and references therein). This is signif- 
icantly motivated by the belief that it could provide an 
answer to the Dark Energy problem and the coincidence 
problem (see e.g. [2-5]). Although it has not been shown 
to be the case, the physics of averaging are still worthy 
of investigation because the parameters of cosmological 
concordance model are quite sensitive to the backreac- 
tion effect [6] such sensitivity is very important in the 
era of precision cosmology. 

One method of evaluating the backreaction effect lies 
within the standard cosmological model. That is, one 
can evaluate the backreaction from perturbations of a 
background spacetime which describe structure forma- 
tion. At second-order, this effect gives rise to corrections 
to the local Hubble flow [6, 7]. This idea was first investi- 
gated in an Einstein-de Sitter model in [3] , and followed 
up in more detail in [8-10]. This study was extended 
to include the case of a cosmological constant in [6, 11- 
13]. However, there appears to be some discrepancies 
between these results: While [8-10] found an important 
effect from backreaction, [6, 11-13] found much smaller 
changes to the value of Hq . Therefore our aim in this 
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paper is to reconcile these results, and present them in a 
unified framework. 

We will make use of the averaging formalism developed 
in [6, 14, 15], to estimate the corrections to averaged lo- 
cal Hubble flow by the small scale inhomogeneities in 
the matter distribution. Already these studies [6, 8-12], 
have used this formalism to evaluate the corrections to 
averaged Hubble rate, averaged deceleration parameter 
and the equation of state of matter fluid, to be expected 
from the small scale inhomogeneities in the matter distri- 
bution, but in each case the definition of averaged Hubble 
rate used is different. For instance, the authors of [8-10] 
, considered a definition of Hubble rate which is very dif- 
ferent from the one considered in [11, 12] and also in [6]. 
The authors in each case considered different slicing of 
the averaging domain, and/or different approaches to the 
cosmological perturbation theory. Specifically in [11, 12] 
, the authors defined the averaged Hubble flow in the 
longitudinal gauge by following the expansion of the co- 
ordinate grids adapted to the gauge; this is an expansion 
rate which is associated with the gravitational potential 
since the magnetic part of the Weyl tensor vanishes in 
this gauge [6] (see also the Appendix of this paper for 
details on the magnetic part of the Weyl tensor). This 
formalism was later used in [12] to study the effect of 
backreaction from averaging in various gauges. 

On the other hand, [8-10] looked at the expansion of 
the matter fluid in the comoving synchronous gauge [8- 
10], while [6] studied the same thing, but in longitudinal 
gauge. Kolb et. al [3] studied the relationship between 
the averaged expansion of matter fluid in the comoving 
synchronous gauge and the longitudinal gauge. The re- 
sults and claims from these papers appear to differ base 
on the approach used. It is against this backdrop that 
we propose to study quantitatively the difference between 
different definitions of Hubble rate. This study will en- 
able us to clarify this issue and at the same time evaluate 
precisely the corrections to the concordance model as a 
result of the backreaction effect. We will also discuss 
how a consistent second order treatment of the backre- 
action effect in cosmological perturbation theory changes 
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the value of the Hubble rate. The analysis of the intrin- 
sic variance which could affect the measurement of the 
Hubble rate today, H , will be thoroughly investigated. 

This paper is organized as follows: In Sections II and 
III, we briefly recall the averaging formalism to be used, 
and discuss the definitions of the averaged Hubble rate 
in two different hypersurfaces of interest. In this section, 
we will also use the gauge invariant formalism developed 
in [16, 17] to calculate the averaged Hubble rate defined 
in the fluid frame. To the best of our knowledge, it is 
the first practical calculation making use of of this gauge 
invariant formalism to study backreaction effect. The 
discussion of our results is presented in Section IV and 
a fitting formula for the variance of the Hubble rate will 
be given here. We also show that the two classes of def- 
initions can be clearly distinguished. A brief comment 
on their relevance is also given. Finally, in Section V, 
we draw some conclusions and discuss future works. The 
Appendix presents the detailed expressions of the various 
Hubble rates considered in this paper, at second oder in 
cosmological perturbation theory. 

Throughout this paper, we will suppose that gravita- 
tion is well described by general relativity on all scales 
and that the cosmic matter fluid can be considered as a 
perfect fluid. Moreover, Latin letters of the beginning of 
the alphabet (a, b, c, h) will denote spacetime indices, 
and Latin letters in the middle the alphabet k, ...) 
will denote spatial indices. 



II. EQUATIONS OF MOTION 

Buchert's averaging formalism [14, 18] (a similar aver- 
aging formalism was presented earlier in [19]) and its gen- 
eralization to arbitrary coordinate systems [12, 13, 15, 20] 
rely on Einstein equations written in the Arnowitt-Deser- 
Misner form. Within this formalism, one considers a 
set of observers defined at each point of the spacetime 
manifold, and characterized by a unit 4-velocity field, 
n a , that is everywhere timelike and future directed, i.e. 
n a n a = — 1, with zero vorticity. This 4-velocity field in- 
duces a natural foliation of spacetime by a continuous 
set of space- like hypersurfaces locally orthogonal to n a . 
The projection tensor field onto these hypersurfaces is 
defined as h a b = g a b + n a nb- The line element can then 
be written with respect to this foliation: 

ds 2 = -{N 2 - N,N l )dt 2 + 2N,dtdx l + h lj dx l dx j , (1) 

where we have introduced respectively the lapse function 
N(x a ) and the shift 3-vector N l (x a ) . The components 
of the 4-velocity of the fluid comoving with the coordi- 
nate grids is given in relation with the lapse and shift 
functions as n a = ±(l,-N l ), n a = N(-l, 0,0,0) .It is 
orthogonal to the hypersurface h a b ■ The intrinsic curva- 
ture of the hypersurfaces is given by 1Z = h ab lZ a b, where 
lZ a b is the 3-Ricci curvature of the hypersurfaces and 
the extrinsic curvature (or second fundamental form): 
K ab = -h%hfn C ;d- 



Here we will consider only the Hamiltonian constraint 
and the evolution equation for the metric of the spacelike 
hypersurface. (for the complete set of ADM decomposed 
Einstein equations see [15]). 

(d t -Cv t )h ij = -2NK ij , (2) 
R + K 2 -K l3 K ij = Wire , (3) 

where e = n a n b T a b and T a b is the energy momentum 
tensor defined to include the cosmological constant as 
Tab = (p + p)u a u b + (p + k/{%irG))g ab . u a is time-like 
4-velocity for the matter field normalized to u a u a = — 1, 
it is related to n a through 

u a = -y(n a + v a ) , where 7 = 1 . (4) 

VI - v a v a 

The vector v a is spacelike and it is orthogonal to n a 
(v a n a = 0). 

The non-local, free gravitational field is described by 
the Weyl tensor. Given a timelike vector this is split into 
electric and magnetic parts. For example, with respect 
to n a these are 

E a f = C acbd n c n d and ffW = *C acbd n c n d , (5) 

where C a bcd is the Weyl tensor and *C a bcd is its dual. 
Analogous definitions exist for the vector field, u a . This 
means that observers in the frame of the fluid and ob- 
servers in the coordinate frame observe this electric- 
magnetic split differently (see [21] for the transformation 
relations between the two), analogously to boosted ob- 
servers measuring different electric and magnetic parts of 
the electromagnetic field. In particular, in certain grav- 
itational fields there may exist a special frame whereby 
one of these two components vanishes. For example, in 
so-called silent universes which are not conformally flat, 
there exists a preferred frame in which the magnetic part 
of the Weyl tensor is zero - such a frame may be consid- 
ered the rest-frame of the gravitational field. In space- 
times where this is possible, it is unique as follows from 
the transformation laws in [21], and there exist (at least) 
two physical, well motivated, frames: the rest-frame of 
the fluid, and the rest-frame of the non-local gravita- 
tional field. 

So far we have defined two different 4- velocities, which 
according to standard 1 + 3 decomposition of a covariant 
derivative of 4- vector, will imply defining respectively two 
expansion rates. 



II. 1. Decomposition of velocities 

The covariant derivatives of the two 4— velocities, u a 
and n a , as well as the spacelike relative velocity v a , may 
be invariantly decomposed with respect to the coordinate 
frame, n a , (this corrects the decomposition presented in 
[15]; however the expression for the Hubble rate is not 
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affected) : 



V a n b = -n a h b + -£h ab + S afc , 

V a u b = -jv c (j 2 v c + h c ) n a n b 

1 



(6) 



-7 ^7 v c V a v c + -£,v a + S ac w c J n b 
+-fn a (-f 2 v c v c v b + n {b) + vq)) 
+ ^9h ab + (T ab + uj ab , (7) 

V a v b = -h c v c n a n b - n a v {b) + Qc«a + Z ac v c ^J n b 

+ ^ Kh ab + Pab + W ab , 



(8) 



where: 



£ = h aD \7 a n b , S ab = h c a h a b V (c n d) - -£h ab , 
9 = h ab V a u b , a ab = h c a hfV (c u d) - ^6h ab 



u ab = h c a hy [c u d] , k = h ab V a v b 



Pab = K h bJ(cVd) - ^nh ab , W ab = h c a hffi [c v d] . 

Where we have used the notation A a ... b — n c S7 c A a ... b , 
the angle brackets denote symmetric, trace free, and pro- 
jected with respect to n„. V denotes the spatially pro- 
jected covariant derivative. 

Here £ and 9 are the expansion rates, while k is the 
divergence of the 3-velocity v a ; T, ab , a ab and (3 ab are the 
shear, while Lo ab and W ab are the vorticity in the respec- 
tive definitions. Every quantity defined here has a natu- 
ral interpretation in terms of observers comoving with the 
fundamental 4— velocity n a . Provided these definitions 
are unique and consistent, all related quantities have a 
direct physical or geometric meaning with respect to the 
fundamental 4— velocity n a . Any difference between such 
4- velocities will be of 0(e) in perturbed FLRW case and 
will disappear in the FLRW limit [21] . Note that: 

9 - 7£ = 7« (l + \^ v2 ) + 7 3 « a «Vafc, (9) 

in FLRW limit, 7 and v 2 are of the order 0(e 2 ), hence 
9 ~ £. The decomposition of the matter 4- velocity, u a , 
is quite unusual, since it is with respect to n a . One can 
calculate directly the normal acceleration, vorticity and 
shear and so on; for us the intrinsic expansion rate is 
important: 



6 = V a u a 

= 9 + "fV a (j 2 V a + h a ) 



(10) 



The following relations between expansion rates will 
be used later: 



- 1~ 1 {< j ij + a Bij) ■ 



(11) 

(12) 



where 9b = — 7« — 7 3 -B and for the shear: 

&Bij = -jfiij - 7 3 ( B (ij) - ^ Bh ij)- Trie tensor B ab is 

defined as 

B a b = 7^n{v a n b + v a v b ) + f3 ca v c n b 

+/3 ca v c v b + W ca v c n b + W ca v c v b , (13) 
and its trace is given by B = \nv 2 + (3 ab v a v h . 



III. AVERAGED HUBBLE RATES 

In general, the average of a scalar quantity S(t, x) may 
be defined as: 



(S(t,x)) v 



Jd 3 xJS(t,x) 
Jd?x~J : 



(14) 



where J = \fh is the square root of the determinant of the 
metric on the hypersurface orthogonal to n a . The time 
derivative of Eq. (14) leads to a commutation relation 
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[Or, (-) v ]S(t, x l ) = (N£S) V - (N0 V (S) v 



(15) 



as is usual in the averaging context. 

There are different definitions of the averaged Hubble 
parameter H-p in the literature, and we would like to 
be able to compare them in the context of the standard 
model, up to second-order in cosmological perturbation 
theory. We shall employ the longitudinal gauge below 
in order to calculate averages in the concordance model, 
which fixes our coordinate frame n a . In the longitudinal 
gauge the magnetic part of the Weyl tensor vanishes, and 
the electric part is a pure potential field in the absence of 
anisotropic stress [6] (see also Appendix A), making this 
the rest-frame of the gravitational field, or Newtonian 
frame. In this sense, both n a and u a are physically well 
defined reference frames. 

There arc different local expansion rates: 

• £: the expansion of the family of coordinate ob- 
servers. In the longitudinal gauge we employ below, 
this is the rest-frame of the gravitational field. 

• 9: The expansion of the fluid, as observed in the 
fluid rest-frame. 

• 9: The expansion of the fluid, as observed in the 
gravitational rest-frame. 

When performing averaging, there are two spatial hyper- 
surfaces of interest: 

• ( ) v : Averaging in the gravitational frame. 

• ( )jr: Averaging in the rest frame of the fluid. 

Finally, when averaging expansion rates associated with 
the gravitational field, there is the issue of the time coor- 
dinate to use: we can associate the time coordinate t with 
the proper time of the 'averaged observers', which, when 
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using n a requires an extra factor of N in the expansion 
rate [6]. 

Definitions based on £ 

As argued in [11, 17], one can consider the evolution 
of the metric of the hypersurface: 



dthij = -NhijZ + 2JV£ tf + DiNj + D N t (16) 

and also assume that the dimensionless domain scale fac- 
tor can be defined as a-p — ( yp(o) ) where Vd is the 
volume of the domain. It is easy to show from equation 
(16), that 

3H V = = J- f (m + D k N k ) Vhd 3 x 

Vt> Vt> Jx> 



(N£ + D k N k ) 



v 



(17) 



This definition describes the average expansion of the 
coordinate grid and says nothing directly about the mat- 
ter field. It has been used in the recent literature for 
calculations in the longitudinal gauge [11, 12], in which 
case it can be interpreted as the expansion rate of the 
gravitational rest frame. We will find that this definition 
exhibits some interesting features, such as weak scale de- 
pendence of backreaction effects. 

With this definition, the averaged Hamiltonian con- 
straint ( 3) becomes: 

ml = 16ttG <iV V (p + v 2 p)) v + 2 (N 2 A) V 

-(N 2 1l) v -Qv+Tv , (18) 

Vv = \ (NOv (D k N k ) v + 2 - (D k N k ) v , 

where Qd is the usual backreaction term and Vt> is an 
additional backreaction term which arises because of the 
inclusion of the shift parameter N k . This definition was 
used in [11, 12]. Most studies have tended to include the 
lapse function in the definition. However, this choice is 
arbitrary (extensive discussion on this was given in [6].). 
One can also choose to define the Hubble factor without 
the lapse function as 3H?> = (£,) v and the corresponding 
averaged Hamiltonian constraint becomes 



ml = 16ttG < 7 2 ( P + v 2 p)) v + 2 (A)^ - (K) v 



-Q 



(19) 



Definitions based on 9 

Assuming all types of matter follow the same 4- 
velocity, the local expansion of the matter is given by 
0. If we average this on spatial surfaces orthogonal to 
u a , we have 3Hx> = (&)jr- This definition is equivalent 
to that studied in [7, 9, 10], and is the same as the ex- 
pansion of the coordinates if we choose the synchronous 
gauge. The equations in that case are well known and 
presented in [1]. This choice might be most natural for 



supernova observations, where the observer here on earth 
is assumed comoving with the source. 
Definitions based on 9 

A final definition of the expansion we consider is given 
by 9: the derivative of the matter observers worldline 
projected into the rest-space of the gravitational frame. 
This was introduced in [6, 15] as a way of recognising 
the fact the rest-frame of the matter before and after 
averaging are not the same. Hence, a useful definition 
of the average Hubble factor is: 3H?> = {N9) v .This will 
lead to the following averaged Friedmann's equation: 

ml = 167TG ((j 4 N*p) v + ( 7 2 ( 7 2 - 1)N 2 P ) V ) 

+2A {N 2 1 2 ) v - (7 2 N 2 K) V -Qv+£v, (20) 

= ! (<(OT) 9 >„ - (N9) 2 V ) -2(N 2 a 2 ) v , 

Cv = 2 (N 2 a 2 B ) v - I ((N9 B ) 2 ) V | (N 2 99 B ) V . 

where a 2 = a^aj /2 and o 2 B = <r B j<J Bi /2 + crijcr B had be- 
ing adopted for simplification purposes. In the same vein, 
we can also consider a definition of average Hubble fac- 
tor without scaling with lapse function as 3Hz> = {9) v , 
in this case the averaged Friedmann's equation becomes 



mi 



16ttG « 7 V) P + ( 7 2 (7 2 - 1)P) V ) + 2A (r% 



- {i 2 n) v - Q v + C v , (21) 

C v ^2(a 2 B ) v - 2 -((9 B ) 2 ) v -^{99 B ) v . 

Notice that the Friedmann part of the Buchert equations 
averaged on the comoving hypersurface may be recovered 
from the last two definitions in the limit where Vi — > 0, 
7 — > 1 and (6) v — > [1]. The equations above [18- 

21] contain the standard backreaction term Qx> and ad- 
ditional backreaction term Cx>- The additional Cx> term 
exist because of the non- vanishing peculiar velocity, v a . 
In almost FLRW metric, its contribution will be subdom- 
inant 0(e 2 ). 



III.l. Spatial averaging of a perturbed FLRW 
model 

The equations derived in Sec. II are not closed, but 
physical information can be extracted from them if we 
suppose that the Universe is well described by a per- 
turbed FLRW background. We shall consider pertur- 
bations in the longitudinal (Poisson) gauge, where the 
metric may be written as 

ds 2 = - (1 + 2$ + $ 2 ) dt 2 + a 2 (1 - 2$ - * 2 ) 5 ij dx i dx j . 

(22) 

Here, the coordinates are chosen to coincide with the n a 
frame such that n a — —Nd a t, where the lapse function is 
N = (l + $ + i$ 2 - j^ 2 )- We have used the trace-free 
part of the momentum constraint to set: *i = <f>i = <f> 
(that is, there is no anisotropic stress at first-order). It 



■5 



was shown in [22-25] that the vector and tensor modes in- 
duced by scalars are subdominant at second order, hence 
this line element is sufficiently accurate for our purposes. 
The peculiar velocity v l can be expanded to second order 
and is given by: 

Vi = ^^( 2v i + v z)- 

As usual, the background Friedmann's equation and the 
deceleration parameter for the pure dust and positive cos- 
mological constant universe are given by: 

H(z) 2 = Hi [O (l + zf + 1 - fi ] , 

. . la „ 1 + z dH „ 3 . . 



respectively, where 



Q m (z) 



n {i + z) 3 



[n (i + z) 3 + i-n ] 



1/2 



and the first order peculiar velocity given in terms $ 
reads: 



V W = 



3aH 2 Q r 



+ (23) 



For details about the solution to the first and the sec- 
ond order equations used in this work, see Appendix A. 

In this framework, the average quantities on the hy- 
persurface orthogonal to n a can easily be expanded to 
second order in perturbation theory, so that one would 
rather evaluate Euclidean integrals instead of a Riemann 
integrals [3]: 

(S) v = S + (S 1 ) + (S 2 ) + ^- (SiJi)--^ <Si) (JO . (24) 

where J and J x respectively stand for the background 
and the first order piece of the square root of the deter- 
minant of the 3— metric, Vh, while So, Si and S 2 are the 
background, first order and the second order component 
of any perturbed scalar on the hypersurface orthogonal 
to n a . Note the important terms of the form ( ) ( ) which 
appear due to the Riemann average - such terms do not 
appear if we average perturbations on the background 
only ( as in [26] for example). 



III. 1.1. Frame switching 

In order to perform spatial averages on the hypersur- 
face comoving with the matter fluid, i.e. on the hypersur- 
face orthogonal to u a , while using the coordinate system 
of the longitudinal gauge presented in Eq. (22) we em- 
ploy the technique developed in [17]. This will allow us 
to perform an average in a frame which is tilted with 
respect to the coordinates. We do this because the lon- 
gitudinal gauge is well defined at second-order, and the 
solutions up to second-order are known in the case where 
the cosmological constant is non-zero [27]. 



Before applying the formalism of [17] to the particular 
case of interest here, we summarise it and generalise it 
for our purposes. When defining the average of a space- 
time scalar there is considerable freedom in the definition, 
and this freedom can be used to switch from an average 
defined in one frame to that in another ([17] used it to 
define gauge- invariant averages). Consider defining the 
average of a quantity using a spacetime window function 



(S) 



f M # X y/=gSW a (x) 



j Mi d±x^—gWn(x) ' 
where a a suitable window function might be: 



(25) 



W n {x) = S(A(x) - A )H(r - B{x)) . (26) 

In this definition of the window function, H is the Heav- 
iside step function and B(x) is a positive function of the 
coordinates with space-like gradient, V a B, and A is a 
suitable scalar field with time-like gradient, W a A, such 
that it takes on a constant value Aq on the hyper-surface 
of interest. The scalar field A then defines the folia- 
tion of spacetime for averaging. The range of integration 
across the hyper-surface is specified by inserting a step- 
like definition of the spatial boundary using the function 
B(x), which is then bounded by a constant positive value 
r > 0. 

It was argued in [16, 17] that one can consistently in- 
tegrate out the coordinate time to define an average of 
the scalar field S on the hypersurface of constant A by 
performing a suitable change of coordinates that trans- 
forms the integration variable from This can be 
achieved by defining t = f(t,x), where the function / 
is chosen to ensure that the scalar field S transforms 
as S(f(t,x),x) — S(t,x). By the use of the Jacobian 
factor dt/dt, the 3- metric is also transformed from hij 
into another metric hij. The function / ensures that 
the scalar field A(x, t) is homogeneous in the tilde frame: 
A(f(t,x),x) = A(t,x) = A(°\t) (see [17] for details). 
Inserting this into Eq. (25), one finds: 



It, A d z xJ S(to,x) 
L d 3 xJ 



(27) 



where the tilde quantities are evaluated in the new coor- 
dinate system. According to [17], this represents a gauge 
invariant prescription for the average of a scalar object 
S on the hypersurface Sa °f constant A = Aq. 

In cosmological perturbation theory, the square root of 
the determinant of the metric g~ a b and the scalar field S 
can be expanded to second order in perturbation theory 
to give the average of a scalar field in the new coordinate 
system: 

(S) Ao = S + (Si) + (S 2 ) + j ( J1S1) - j (Si) ( Ji 

(28) 

where Jo = y/—go an( l Ji = V~9i are the background 
and the first order piece of the square root of the determi- 
nant of the metric g a b respectively. The metric g is still 
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a 4-dimcnsional metric. By making a gauge transforma- 
tion [28] back to the original coordinates, we obtain: 



(S) Ao =S q + (Si) + (S 2 ) + -1 (JiSi) - 4" {Si) (Ji) 



5o_ 



(Ai) + — (Ji^i) + (A 2 ) 
Jo 



1 



+: 



1 



So o A Sv dt (In Jo) Sq ■ 

P ~ ~aF~ — P — 



A 1 S 1 ) + (S 1 A 1/ 

{A) 



A d t (In Jo) 



A 

A d t (lnJ ) 



(SiAi) + (Ai) (Ji) 

Jo^o 



Ao 

S*oA , 9 t (lnJ )5 



(Si) (^i) 



Al 



Al 



(Ai) 



(29) 



Here Jo = \f—g~a and Ji = \f—g~i are the background and 
the first order piece of the square root of the determinant 
of the four dimensional metric g respectively. This is the 
major difference between this averaging prescription and 
the conventional one defined in equation (24). Once the 
scalar variable A is chosen to specify the averaging hy- 
persurface, the above averaging prescription can easily 
be applied. Eq. (29) was first derived in [17], but the au- 
thors set the spatial average of a first order scalar quan- 
tity (Si) to zero (see Eq. (3.10) in [17]) before performing 
the gauge transformation, thereby neglecting the terms 

of the form (Si) (A\), (Si) (^Aij, etc, which are non-zero 

and are explicitly scale dependent at second order [6]. 
We have inserted them as they play an important role in 
the average of the Hubble rate. 

To fix the definition of A in terms of the quantities of 
the perturbed FLRW background and at the same time 
fix the foliation of interest, we employ the technique used 
in [3]. This involves relating the scalar field A to the time, 
r, measured by the observers with 4- velocity u a comoving 
with the fluid: u°do + u l di = d T . The scalar field A can be 
expanded to second order in perturbation theory, subject 
to the condition A(t, x) = Ao(t) + Ai (t, x) + A 2 (t, x) = r 
[16] to give (using u a V a r =1): 



(1 - $ - — <1>2 
1 2 



-$ 2 +v^v lk )d t A(t,x) 



(30) 



+ -{v\+v\)d l A(t,x) = l 



We can now calculate the higher order A in terms of $i 



and $ 2 of the perturbed FLRW background. This gives: 
A (t)=t, (31) 



$idt, 



Ai(t,x) = f 
Jo 

A 2 (t,x) = ^J^ <f> 2 dt-^J^ $ldt- J v\vi k dt 



(32) 



T v{diAidt 



(33) 



The average Hubble factor calculated using this prescrip- 
tion is given in the Appendix B. 



III. 2. The ensemble average and the variance 

With the tools developed in Section. III.l, we have 
performed a consistent second order perturbative expan- 
sion of the Riemann average defined in Sec. II. 1 to obtain 
a corresponding Euclidean average. Given a specific re- 
alisation of a cosmology, we could now calculate spatial 
averages directly. Alternatively, we can calculate the en- 
semble average of a given spatial average which will tell 
us the expectation values of spatially averaged quanti- 
ties. The ensemble-variance tell us how much we expect 
that to vary from one domain to another. 

The ensemble average of a spatial average may be writ- 
ten as: 



d 3 xW(x/R v )X(x) 



(34) 



where the overbar denotes an ensemble average. We have 
specified the domain though the window function W. 
The Euclidean volume of the spatial domain of averag- 
ing T> is then given by: V — J d 3 xW(x/Rx>) which in 
the case of a Gaussian window function which we mostly 
employ is V = 4nR^ y 2 W(y)dy = (2n) 3 / 2 R 3 , for any 
Rv- The inverse Fourier transform of this window func- 
tion reads: W(kR v ) = ± J d 3 xW(x/R v ) e - lk x . The 
Fourier and the inverse Fourier transforms of any scalar 
quantity $ are given as 



$(x) = 



*(fc) = 



1 



(2tt) 3 / 2 

1 

(2tt) 3 / 2 



d 3 k<S>(k)e 
d 3 x$(x)e~ 



ik-x 



(35) 
(36) 



For statist ically homogenous Gaussian random variables, 
we have: <I>(fc) = 0, and the power spectrum of $ is 
defined by 



2tt 2 



$(*)$(*') = —V^(k)S(k + k'). (37) 

Assuming scale- invariant initial conditions from inflation, 
this is given by 



(38) 
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where T(k) is the normalised transfer function, A 2 ^ is 
the primordial power of the curvature perturbation [29] , 
with A^ w 2.41 x 10~ 9 at a scale k C MB = 0.002Mpc _1 
(for the definitions of g{z) and g^, please see equation 
A5 in Appendix A). 

It is not difficult to notice from the equations displayed 
in the appendix that most of the terms we are dealing 
with are scalars which schematically appear in the form 
d m <&(x)d n <&(x) where m and n represent the number of 
derivatives (not indices) , such that m + n is even so that 
there are no free indices. (For example, di$>d 2 d l & has 
m = 1 and n = 3.) Then from the results of [6], the 
ensemble average of these kind of terms, if a Gaussian 
window function is assumed, is given by: 



(d m <f>(x)d n $(x)) 



(_l)(™+3n)/2 

2^2 



/ 



dk fc m+ "- 1 fc 3 P$(fc). 

(39) 

Using $ = g(t)$ (x), g(t) being the growth suppres- 
sion factor and &o(x) the spatial dependent initial con- 
dition (see the Appendix), the terms that appear with a 
time derivative of the gravitational potential can be re- 
written to pull out the time component before evaluating 
the ensemble average: 

*(t,x) = -{l + z)H^^*(t,x). (40) 
az 

For the details of the calculation of the ensemble aver- 
age of the inverse laplacian appearing the second order 
Bardeen potential refer to [6]. 

The ensemble variance in the Hubble factor is given by 



Var[ff*] =H%-H X , 



(41) 



where Hx can be any definition of averaged expansion 
rate we are studying. With this definition, it is easy to 
see that pure second order contributions drop out of the 
variance, so that only terms that are quadratic in first 
order quantities remain. 



IV. RESULTS AND DISCUSSION 

We shall now investigate the expectation values of the 
different average Hubble rates, along with their variances. 
For this we will consider an Einstein-de Sitter model, and 
a standard concordance model. We shall use length scales 
intrinsic to the model as reference points for averaging: 
these scales are the equality scale, A;" 1 , and the Hubble 



k cq « 7.46 x I0- 2 n h 2 } 

fc H ~ Mpc \ 

3000 P ' 



(42) 



where fit and Slo are the baryon and total matter contri- 
butions today and H a = 100 /ikms _1 Mpc~ 1 . We shall 
also use two models for comparison. The first is an 
Einstcin-de Sitter model with h = 0.7 and 5% baryon 



fraction (WMAP [29] whose estimate as an energy den- 
sity is given as 6 w 0.046). This has fc^ 1 ~ 27.9Mpc. 
The other model we shall use is the concordance model 
with Sl n = 0.26, h = 0.7, /baryon = 0.175 (this is the 
WMAP best fit [30]). The key length scales in these 
model are fc^ 1 ~ 107.2Mpc and the Hubble scale k^ 1 ~ 
4.3Gpc. 

To calculate the integrals we shall use the transfer func- 
tions presented in [31]. All lengths scales shown are in 
Mpc unless otherwise stated. Because some of the in- 
tegrals have a logarithmic IR divergence, all fc-integrals 
have an IR cut-off set at ten times the Hubble scale, 
it did not appear explicitly in any of our calculations. 
Moreover, the position of the IR cut-off does not affect 
the result, that is one can set the IR cut-off within this 
range 10 — 10 9 and the results shown here will remain un- 
changed. Since the divergence is logarithmic, the result 
depends very mildly on where the cut-off is set. 

We show the ensemble averages of some of the second- 
order terms which appear in the Hubble rates in Fig. 1. 
Note that we also show the result of (9 2 $<9 2 $) for com- 
parison. 



IV. 1. Comparison between the different definitions 

We can now turn to estimating and comparing the 
Hubble rates as well as their intrinsic variances as defined 
above consistently up to second order in perturbation 
theory. When determining the ensemble average of the 
Hubble rate, we shall consider two alternatives: a kine- 
matical ensemble average given by H-u, and a dynamical 
one, which arises from taking the ensemble average of the 

Friedman equation: \J iJ 2 ,. We shall find that the differ- 
ence between these two is large because the variance is 
large. 

Fig. 2 presents the evolution of the averaged Hubble 
rate as functions of redshift for different definitions of 
Hubble rate in a ACDM and an EdS scenarios. Fig. 3 
depicts the values of the same Hubble rates at z = as 
functions of the averaging scale Rd, and Fig. 4 shows the 
scaling of their variances with the averaging scale Rd. 

It is clear that the two types of Hubble rates defined 
here, i.e. those of the gravitational frame, and the ones 
defined in terms of the physical matter flow can be dis- 
tinguished as far as the magnitude of their mean and 
variance are concerned. 

First, the ones defined through the local expansion 
of the observers' worldlines, and (N^) v present a 

very small correction to the FLRW background Hubble 
rate, which is of the order 10~ 5 for ACDM and 10" 4 
for an EdS scenario at 20 ft, -1 Mpc . Such a small effect 
was also reported from a study using numerical simula- 
tion [32]. Moreover, they appear to be scale independent 
when compared with the definition based on 6. The scale 
dependence is determined by terms with two angle: brack- 
ets (<&) ($). For the definitions based on £, for example, 
the variance in (£) v or (N£) v is scale dependent but 
its scale dependence is suppressed when compared with 
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FIG. 1. The amplitude of the ensemble averages of various 
second-order terms which appear in the Hubble rate. 
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FIG. 2. Fractional change to the background Hubble rate 
as a function of redshift for the different definitions of av- 
eraged Hubble rates under study. Here we have averaged 
at the equality scale. Here, Ah = (-Hp — Ho) /Ho, and 

$h = {JhJ-H )/H . 



in £ definition is ($) ($) which is of the order of 10 10 
at 20 /i _1 Mpc while in 6 definition the dominant term is 
(<9 2 $) (<9 2 $) which is of the order of 10~ 4 at 20 fi^Mpc. 

Second, the Hubble rates defined through the local ex- 
pansion of the matter worldlines systematically present 
corrections to the background Hubble rate, which is two 
orders of magnitude bigger than the previous ones, and 
are indistinguishable from each other, except when the 
averaging scale is much larger than the equality scale. 
It is interesting to note that both the values of these 
averaged Hubble rates and their variances are indistin- 
guishable up to scales of averaging of order >100 Mpc, 
after which they start to differ. This scale have been 
interpreted in a previous work [6] as naturally defining 
the scale of statistical homogeneity of the universe (note 
that this is the case even for EdS; so it is not simply 
the equality scale) . Around the same scale the expansion 
rate of the gravitational frame becomes comparable with 
the others because the peculiar velocity tends to zero. 

Finally, let us note that the results are consistent, for 
a pure CDM Universe, with those found on small scales 
in [9, 10]. This can be seen on Fig. 5. 

This analysis shows that the averaged Hubble rates 
defined through the expansion of the Newtonian-like or 
gravitational frame, as in [11, 12], is not a good tracer 
of the expansion of the cosmic fluid, except beyond the 
homogeneity scale. The fluid frame is more relevant for 
local measurements since it is attached to the matter 
component of the Universe. The 'gravitational frame', as 
we have referred to it here, seems useful on much larger 
scales, which is the situation in [11, 12] in which it was 
first evaluated - their domain was the Hubble scale. 
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FIG. 3. Fractional change to the background Hubble rate as a 
function of the averaging scale for the different definitions of 
averaged Hubble rates under study. In both models, (£),£, and 
{NQv almost coincide and are indistinguishable in the figure. 
Notice the turn-down in (Q)j- definition in Einstein de Sit- 
ter Universe, the physical meaning of this is not immediately 
obvious 



IV. 2. Fluctuations in the measurement of H { 



o 



We would like to finish this paper by addressing the 
following questions: 

• What is the physical relevance of the averaged Hub- 
ble rate and its variance? 

• Can there be any signature of backreaction in the 
observations leading to the measurement of Hq? 

First, let us note that on sufficiently small scales, such 
as scales smaller than ~ 100 Mpc, which are the stan- 
dard scales at which the Hubble rate is evaluated, and in 
a statistically homogeneous and isotropic Universe, spa- 
tial averages are expected to be a good approximation of 
what happens along the past lightcone on which observa- 
tions are made. Along the past lightcone the monopole 
contribution to the Hubble rate, which is the one that 
remains once a full sky average has been performed, is 
exactly the covariant quantity 6 = V a u a [33]. 

Hence, our estimate of H-p on a scale Rt> can be inter- 
preted as the average Hubble rate in a patch of the local 
Universe of size Rt> as long as this size remains suffi- 
ciently small compared with the Hubble scale. Moreover 
the variance we calculated is the intrinsic dispersion on 
the measurement of H that comes from the fluctuations 
in the peculiar velocity of the sources and gravitational 
potential. In a concordance cosmology, this dispersion 
appears small, of order 1% at a scale of 100 Mpc, and 
even less on larger scales, as can be seen on Fig. 5. 
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FIG. 4. Variance in the fractional change to the background 
Hubble rate as a function of redshift for the different defini- 
tions of averaged Hubble rates under study. 



This is consistent with previous estimates that were 
based on an estimate of the first order velocity power 
spectrum [7, 34]. It is due to the fact that the pure second 
order terms cancel out consistently at second order in our 
expression of the variance, allowing only contributions of 
squares of first order quantities. As noted before, this is 
a similar effect to that found in [9, 10], where the calcu- 
lations were made in the comoving synchronous gauge, 
for a pure CDM Universe. Our calculation of (6)jf us- 
ing the gauge-invariant approach of [16] corresponds to a 
gauge-invariant version of the average expansion rate in 
the synchronous gauge. 

To quantify the backreaction effect on the variance for 
a large class of cosmological models, we provide a fitting 
formula for the variance of the Hubble rate (defined via 
the flow of matter), Var[H], that is accurate to a few 
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R nT (h" r Mpc) 



found that the definitions that involve the flow of the 
dust matter component are consistent with each other at 
second order in cosmological perturbation theory, but dif- 
fer significantly on small scales from the definition based 
one the expansion of the coordinate grids. In particular, 
we noticed the following features of the averaged Hubble 
rate: 

• On small scales all definitions which involve the 
matter flow agree, and give a small sub-percent 
change to the background Hubble rate. 

• The variance in the average of the expansion rate of 
the gravitational frame are very small and appear 
to exhibit weak scale dependence when compared 
with the definition involving the matter flow. 



FIG. 5. Dynamical Hubble rate {O)^ today, plus/minus the 
variance as a function of the averaging scale (normalized to 
the background Hubble rate), where we have used a top-hat 
window function to define our domain for comparison with [9, 
10]. The blue curve represents the variance calculated in [9, 
10]; differences for small domains are a consequence of the 
different transfer functions used here. 

percents across the scales of interest: 

log w (Var[H]) = -43.61 + 46.0D™ 293 - 0.79 69/ b ^) 
+A(log 10 R) a + 7 exp(-/3(log 10 R) 2 ) 

where 

A = 10.32 - 9.084^ 0469 - 3.61 l// b 00497 (44) 
7 = 1.309 - 2.355fi° ra 055 - 1.073/ b 21778 
/3 = -1.805 + 3.260ft™ 279 - 0.7180/ b a665 
a = 1.222 + 0.0334ft 3 „ 635 + 0.0591/ b 3944 . 

This formula gives the variance on the measurement of 
Hq, normalised to the value of Hq for H-p = ■§ (N9) v 
definition (The fitting formula for other definitions exist 
but they are more complicated and less accurate): fl m 
is the CDM density parameter, / b the baryon fraction, 
and R the length characteristic of the survey, i.e. the 
distance to the farthest object (in units of Mpc). Note 
that this fitting formula is valid for a Gaussian window 
function. Top-hat window functions generically lead to a 
slight increase of the variance. 

V. CONCLUSION 

In this work, we have presented the first comparison 
of the different definitions of the averaged Hubble rate 
that can be found in the literature. We did this by cal- 
culating these various definitions consistently to second 
order in cosmological perturbation theory. Also for the 
first time we have calculated the average of the expan- 
sion rate using the formalism of [16], which the authors 
claim is gauge invariant at some limit. We have also 



• On large scales (much larger than the equality 
scale) all definitions become scale invariant once 
their ensemble average is evaluated. 

• The hypersurface used in averaging is not really im- 
portant when computing the variance for perturbed 
FLRW model, the differences only show up on large 
scales and it is only noticeable in Einstein de Sitter 
models. 

• Including N in the definition of the averaged expan- 
sion leaves a residual effect on large scales, and it 
tends to reduce the backreaction effect. The inclu- 
sion of the lapse function is made compulsory if one 
wants to keep the coordinate time t as the proper 
time in the averaged model (there is a discussion 
on this issue in [6]). But this is only one possible 
choice, since no-one knows how to explicitly con- 
struct the average model in this setting where only 
the scalars are averaged. 

We have also derived the dispersion affecting the Hub- 
ble rate and arising from the peculiar velocities of the 
matter flow. We found an effect consistent with previ- 
ous estimates from backreaction in the literature [9, 10], 
and our results are consistent with effects evaluated pre- 
viously [6, 7]. 

We close with a comment on the origin of the scale de- 
pendence of the various quantities. The scale dependence 
we have found here comes only from 'non-connected' 
terms such as ($) (<9 2< I > ) since the domain size factors 
out of all other terms (for details on how this type of 
terms appear see [6]). Non-connected terms only arise 
when we perform averages in the spacetime itself, which 
many authors on backreaction have stressed is important. 
It is interesting to note that these terms (i.e those ones 
involving laplacian of a gravitational potential, for exam- 
ple (<E>) ((9 2 <I>)) do not appear if we treat perturbations 
as fields propagating on the background, and calculate 
average quantities only with respect to the background 
geometry i.e., if we perform a Euclidean average and not 
a Ricmannian one [26] . 
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Appendix A: Second-order perturbation theory 



The Poisson gauge is particularly elegant for scalar perturbations because with n a defined orthogonal to the spatial 
metric hij , the electric and the magnetic part of the Weyl tensor becomes 



4° = o. 



V a $V,,$ - V a *V b * } (Al) 

(A2) 



In the rest frame n a , then, the gravitational field is silent, and, with ^ = $ is a pure potential field. Hence, n a 
may be considered as the rest-frame of the gravitational field, or the Newtonian-like frame, and so defines natural 
hypersurfaces with which to perform our averages. By contrast, in the frame u a the Weyl tensor has non-zero H a b [6]. 

The Einstein Equation for a single fluid with zero pressure and no anisotropic stress = <&, and $ obeys the 
Bardeen equation 



$" + 3H$' + a 2 A$ = = <£ + 4i?$ + A$ . 



(A3) 



and ' = d/dr), and H = a' /a is the conformal Hubble rate. All first-order quantities can be derived from $. The 
solution to the growing mode of the Bardeen equation may be written as 



$(r],x) = g(r))$ (x) 



(A4) 



where &o(x) is the Bardeen potential today (rj = rjo, z = 0) and g(rf) is the growth suppression factor, which may be 
approximated, in terms of redshift, as [35, 36] 



g( z ) = ^9oo^ m (z) \ n m (z) 4/T - n A (z) + 



1 + Ul m {z) 



1 + 



(A5) 



and is chosen so that g(z = 0) = 1. 

The second-order solutions for Vl^ 2 ) and $( 2 ) are given by [27]. We quote their results directly: 
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* (2) (»7,a0 = B\ (rf) - 2g( V )g m - — (a n] - l)g{r,)g m ) $1 + ( B 2 (rj) - -g(r,)g 



V- 2 (d*$ d t <S> ) 



(A6) 



<f (2) (^) = (b^tj) + Ag 2 {r,) - 2g( V )g m - y(a nl - l)g(v)9m^ H + 



B 2 (v) + o9 2 {ri) efa) + - - -g(v)g 



V- 2 (d^od^) - 3V- 4 d^ (FQodjQo) 



+ B a ( V )V- 2 d^(d l ^ d^o) + B 4 (r,)d l -QodiQo , (A7) 



where Bi(rf) — Hq 2 (fa + 3Qo/2) Bi(rf) with the following definitions 



B^n) = f dvH 2 m.f(v) - i) 2 c(v,v), b 2 (v) 

Jr\m 

rv / 3 



dw 2 (v) 2if(fj)-iy-3 + 3n m (fj) c^fj), (as) 



B 3 (v) 



'In 



dfj e(fj) + - C(r),rj) 



BM 



dfjC(r],fj,) 



and 



C( V ,rj)=g'(rj)a(^g(r ] )n(v)-g(rj)^n( V ) 



(A9) 



(A10) 
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with e{jf) = f 2 (i])/tt m (r)) and 



m - 1 + $Hry < A11 > 



g TO denotes the value of g(r) m )i deep in the matter era before the cosmological constant was important. We also have 
a n \ which denotes any primordial non-Gaussianity present. We set this to unity, representing a single field slow-roll 
inflationary model. For details on how the spatial average of the second order Bardeen Potential may be evaluated 
see [6] 



Appendix B: Hubble rates 



In this appendix, we present the different Hubble rates, consistently at second order. The superscript determines 
the quantity that has been averaged to define the average Hubble rate. 



= H - - 3 ^ ($) + <^$$) - - 
Hlh=H-($)-H{$)-3 ($\ ($) + 2 /$$\ - 3H ($) 2 + - ($ 2 ) 



2 

9 I 



(Bl) 
(B2) 



H%° =H- - 3 (*} ($) + 2 - - (& 2 ) r 



(1+f) «, k _ 2{l + zf 

6 0k 2 m m m 



<9 2 i>) + h (a 2 *) 



+ 



(i + ^) 2 

2(l + z) 2 



-h \i + -z.) (9 fc $a fe $) + -h 2 (i + n m ) (d k $d k <s>) - - (d k $d k i> 

| <J$5 2 $) + ($9 2 $) - ^<9 2 $^> (*) - H (<9 2 $) ($) 



(B3) 



i/|> = H — ^ - H <$) - 3 ^ ($) + 2 (^$$^> - 3ff ($) 2 + -if ($ 2 ) - i 

\2 



(1 + * ) -<ft«*> 2(1 + 2 



i, 'Ml,,, IV 1 

3 



<9 2 $) + J ff<9 2 $) 



2(1 



2 



i + 2^mJ + # 2 (a fc $a fc $) (i + 3a, 

#($<9 2 $) -H(d 2 $>) ($)] . 



2(l + z) 

3n„,H 2 



2 



9fe$9 fc $) +2F (5fc$a fe $ 



$a 2 $) - (<9 2 $) ($) 



H$ = H -(&) -{<Z>)H 



1 + -tl m Hgi 



2(1 + z) 2 

9 ^"^TT) 



<9 2 $) + #(<9 2 $) 



ff($ 2 ) -3ff($) : 



(l+ gi H) + -n m H gi (l-H gi ) 



-2($$) + il±fl(^) 
3 ($)($> + 



-- <$ 2 ) n OT H 



1 - U - \h 9i 



2(1 + ^) 2 
27n 2 H 3 



d k <i>d k <i>) + 2H i + -n m ) (s fc $9 fc $ 



(B4) 



i/ 2 (l + 9fi m )(d fc $d fe $)] 



2(1 + z) 2 



i - ^$a 2 $) + f i - ~# <^ # ($a 2 $) 



1 - \HgA (d 2 $>) (*) - f 1- |iTflrj) ff(c> 2 $) (*) - fi m ff 2 fl/ ($a 2 $) 



(A 2 ) U 2 tt m5 / , (B5) 



where gi = —kr f g(t')dt' . The analytical computation in this paper was performed with the help of xAct [37]. 

I 



[11 T. Buchert, Gen. Rel. Grav. 40, 467 (2008), 


[21 S. Rasanen, 


JCAP 


0402, 003 (2004), arXiv:astro 




arXiv:0707.2153-[gr-qc]|. 


ph/0311257 









13 



E. W. Kol b. S. Matarrese, A. Notari. and A. Riotto 
Phys. RevJ D71, 023524 (2005), |arXiv:hep-ph/0409038[ 
S. Rasanen , JCAP 0611, 003 (2006), |arXiv:astro"- 



ph/0607626 



E. W. Kolb, S . Matarrese. and A. Rio tto. |New J. Phys 
8, 322 (2006), |arXiv:astro-ph/0506534 



C. Clarkson, K. Ananda, and J. Larena, | Phvs. R ev 
083525 (2009) " 



arXiv:0907.3377 -[astro-ph.CO 



Y. Wang, D. N . Spergel. and E. L. Turn er, Astrophys. J 



498, 1 (1998), kXiv:astro-ph/9708014 . 
N. Li and D. J. Schwarz, phys. Rev.|D76, 083011 (2007), 



arXiv:gr-qc/0702043 



D80, 



N. Li and D. J. Schwarz, Phys. Rev. D78, 083531 (2008) 



arXiv:0710.5073 -[astro-pE 



N. 



Li, M. Seikel, and D. J. Schwarz, fcbrtsch. Phys 56 



465 (2008), |arXiv:0801. 3420 --[astro-phj. 



[22 
[23 

[24 
[25 
[26 
[27 
[28 



I. A. Brown, G. Robbers, and J. Behrend, JCAP 0904 



016 (2009). arXiv:0811.4495-fgr-qc 



129 



I. A. Brown, I Behrend, and K. A. Malik, JCAP 0911, 



027 (2009), |arXiv:0903.3264-[gr-qc 
J. Behrend, I. A. Brown, and G 
013 (2008) 



arXiv:0710.4964 



T. Buchert, 



qc/9906015 



Gen. Rel. Grav 



Robbe rs, 
astro- pli] . 



JCAP 



0801. 



32, 105 (2000), ttXiv:gr 



Larena, 



arXiv:0902.3159 
M 



P hys. Rev 

gr-q4 

Marozzi. 



D79, 084006 (2009), 



G~ 



and G. Veneziano, 



Gaspcrini 

1002, 009 (2010), |arXiv:0912.3244'[gr-qc 
M. Gasperini 



JCAP 



qc/0102049 



JCAP 



, G. Marozzi. and G. Ven eziano. 
0903, 011 ( 2009) , larXiv :090 1.1303' [gr-qc]| . 
T. Buchert, |Gen. Rel. Gray] 33, 1381 (2001), ^rXiv:gr- 



[30 

[31 

[32 

[33 
[34 

[35 

[36 



H. R uss. M. H. Soffel, M. Kasai. and G. Borne r. [Phys 
Rev.|D56, 2044 (1997), kXiv:astro-ph/9612218 



T. H.-C. Lu, K. An anda. and C. Clarkson. |Phvs. Rev" 



D77, 043523 (2008), fa,rXiv:0709.1619 -[astro-phi 
T. H.-C. L u. K. Ananda, C. Clarkson, and 
R. Maartens, [JCAP| 0902, 023 (2009), * Brief en- 



try *, a,rXiv:0812.1349-[astro-ph] 



K. N. Ananda, C. C larkson, and D. Wands, Phys. Rev 
D75, 123518 (2007), ^,rXiv:gr-qc/0612013 ~[gr-qc]. 



S. Matarrese, S. M ollerach, and M. Bruni, |Phys.Rev 
D58, 043504 (1998), [a,rXiv:astro-ph/9707278 -[astro-ph] 



D. Baumann, A. Nicolis, L. Senatore, and M. Zaldar- 



riaga(2010), |a.rXiv : arXiv: 1 004 . 2488 - [astro-ph . CO] . 

N. Bartolo, S. Matarrese, and A. Riotto, JCAP 0605, 



010 (2006), arXiv:astro-ph/0512481 . 

M. Bruni. S. Matarres e. S. Mollerach. and S . Sonego. 



Class. Quant 
qc/9609040[ 



Gravj 14, 2585 (1997), |a.rXiv:gr- 



E. Komatsu et al. (WMAP Collaboration), |Astro 



phvs.J.Suppl 



-[astro-ph 



180, 330 (2009), 



arXiv:arXiv:0803.0547 



J. Dunkley et al. (WMAP). |Astrophvs. J. S uppl 
306 (2009), |arXiv:0803.0586 -|astro^phil 
D. J. Eiscnstcin and W. Hu, Astrophys.J 



180. 



(1998), [irXiv :astro-ph /9709 1 1 2 - [astro-ph] 
X. 



Zhao 



and 



496, 605 



arXiv:0912.4750-[astro-ph.Cql 



J. Kristian and R. K. Sachs 
X. Shi and M. S. Turner 



Mathews(2009), 



Apr 



pj 143, 379 (Feb. 1966). 



pj 493, 519 (Jan. 1998), 



arXiv:astro-ph/9707101 . 
O. Lahav, P. B. Lilje, J. R. Primack, and M. J. Rees, 
Mon. Not. Roy. Astron. Soc. 251, 128 (1991). 



S. M. Carroll, W. H . Press, and E. L. Turner, Ann. Rev 



Astron. Astrophys. 30, 499 (1992) 



S. Rasaneii. [JCAP| 1003, 018 (2010), ^rXiv:0912.3370 



[astro- ph.CO] " 
R. Maartens, T. Geb bie, and G. F. R. Ellis, ¥hvs. 
D59, 083506 (1999), [irXiv: astro-ph/9808 1 63f ~ 



[37] D. B rizucla. J. M. M artin-Garcia. a nd G. A. Mena Maru 
gan,|G~en.Rel.Grav.]41, 2415 (2009), |arXiv:0807.0824-fer- 



qc 



Rev 



A UUM universe 




